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Numerical study of the Schwinger-Dyson equation (SDE) for the fermion propagator (FP) to 
obtain dynamically generated chirally asymmetric solution in an arbitrary covariant gauge £ is a 
complicated exercise specially if one employs a sophisticated form of the fermion-boson interaction 
complying with the key features of a gauge field theory. However, constraints of gauge invariance 
can help construct such a solution without having the need to solve the Schwinger-Dyson equation 
for every value of £. In this article, we propose and implement a method to carry out this task in 
quenched quantum electrodynamics in a plane (QED3). We start from an approximate analytical 
form of the solution of the SDE for the FP in the Landau gauge. We consider the cases in which 
the interaction vertex (i) is bare and (ii) is full. We then apply the Landau-Khalatnikov-Fradkin 
transformations (LKFT) on the dynamically generated solution and find analytical results for arbi- 
trary value of £. We also compare our results with exact numerical solutions available for a small 
number of values of £ obtained through a direct analysis of the corresponding SDE. 

PACS numbers: 11.15.Tk,12.20.-m 



I. INTRODUCTION 

Quantum electrodynamics in a plane (QED3) con- 
tinues to attract attention both in the field of super- 
conductivity, e.g., where it has been used in the 
study of high T c super-conductors, as well as in the realm 
of dynamical generation of fundamental fermion masses 
where the numerical findings on the lattice and the re- 
sults obtained by employing Schwinger-Dyson equations 
(SDE), @], are yet to arrive at a final consensus. 

In this paper, we take up the latter question in the light 
of SDE. Owing to its simplicity as compared with its four- 
dimensional counterpart and quantum chromodynamics 
(QCD), the quenched version of QED3 is particularly 
neat to unfold the intricacies of the SDE for the fermion 
propagator (FP). It provides us with an excellent labo- 
ratory to study one of its unwelcome features, namely, 
the lack of gauge invariance of the associated physical 
observables. It hampers a fully reliable predictive power 
of the said equations in the non-perturbative regime of 
interactions. This problem can be traced back to not 
employing, or doing so incorrectly, the gauge identities 
such as Ward-Green- Takahashi identities (WGTI) Q , the 
Nielsen identities (NI) |j] and the Landau-Khalatnikov- 
Fradkin transformations (LKFT) 5]. In this article, we 
address this issue in the light of the LKFT. 

The LKFT of the Green's functions describe the spe- 
cific manner in which these functions transform under a 
variation of gauge. These transformations are nonper- 
turbative in nature, and they are better described in co- 
ordinate space. It is primarily for this reason that some 
earlier works on its implementation in the study of the 
FP were carried out in the coordinate space, In the 
context of gauge technique, the FP in QED4 was found 
to have correct infrared and ultraviolet behaviour on em- 
ploying an ansatz for the longitudinal vertex. Momentum 



space calculations are more tedious, owing to the compli- 
cations induced by Fourier transforms. These difficulties 
are reflected in [8j where non perturbative FP is obtained 
starting from a perturbative one in the Landau gauge in 
QED3 and QED4. 

In perturbative calculations of the Green's functions, 
these transformations are satisfied at every level of ap- 
proximation. However, in the non-perturbative study 
of gauge field theories through SDE in an arbitrary co- 
variant gauge, one can ensure these transformations for 
the FP and the vertex are satisfied only with a correct 
ansatz for the so-called transverse part of the three point 
fermion-boson vertex. In practice it implies constructing 
a very complicated vertex if one wants to obtain correct 
gauge covariance in all momentum regimes. In QED3, 
recently an ansatz has been proposed, [(| which guar- 
antees that the resulting FP satisfies its LKFT to two 
orders and the vertex itself to the first order in their 
respective perturbative expansions. As LKFT are non 
perturbative in nature, we expect them not only to be 
satisfied at every order in perturbation theory but also 
in phenomena which are realized only non perturbatively, 
such as dynamical chiral symmetry breaking (DCSB). 

The continuum studies of DCSB are carried out 
through the SDE for the FP. Gauge dependence of the FP 
can be obtained by solving these equations in different co- 
variant gauges. In practice, one needs to go in small steps 
of the gauge parameter away from the Landau gauge and 
it is prohibitively difficult to be able to compute the re- 
sult for an arbitrarily large value of the gauge parameter 
especially if a sophisticated form of the three point inter- 
action is taken into account, The LKFT provide 
us with a mechanism to achieve this goal [llj . What we 
have to do is to start from the result in the Landau gauge 
and simply perform an LKFT to find the result in any 
other gauge. 
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The paper is organised as follows : In Sect. II we recall 
the solution of the FP in the Landau gauge employing 
the bare vertex. In order to get better physical insight, 
we approximate the solution with a simple analytic para- 
metric form; in Sect. Ill, we perform the LKFT of the 
dynamically generated to get the FP in an arbitrary co- 
variant gauge. In Sect. IV we discuss in details various 
limiting cases of our results. In Sect. V, we go further 
to implement the full vertex. We present the numerical 
results in Sect. VI and conclusions in Sect. VII. 



II. THE FERMION PROPAGATOR FOR THE 
BARE VERTEX 

A starting point for the LKF transformations is knowl- 
edge of the Green's functions, here the FP, in the Landau 
gauge. We begin this study using the rainbow approxi- 
mation (where the vertex is bare) since in this simple case 
the Landau gauge FP is easy to find. The well known 
SDE for the FP in quenched QED3 in the Minkowski 
space is : 



SmHp;0 = S M \p) 



d 3 k 



r» M (k,p)S M (k;Z)rAZ»(q),(l) 



where q = k—p and e is the dimcnsionful electromagnetic 
coupling. The bare fermion and photon propagators, 
respectively, are S%(p) = V^> A fi°(<l) = -9^/q 2 + 
(1 — Ciq^qv/q 4 " and £ is the covariant gauge parame- 
ter. T^j(k,p) is the full fermion-boson vertex which we 
take to be bare, i.e., 7^, for the time being. We shall 
undo this assumption in Sect. V. and implement the 
full vertex. We prefer to write the FP in its most gen- 
eral form as S M (p;0 = F(p;0/(p' - M(p;0). F(p;0 
is referred to as the wavefunction renormalization and 
M(p;£) as the mass function. The subscript M stands 
for the Minkowski space. Eq. Q is a matrix equation. 
On multiplying it by 1 and p 1 and taking the trace, it 
decomposes into a pair of coupled equations involving 
F(p; £) and M(p;£,). After a Wick rotation to Euclidean 
space, and performing angular integration, it is a com- 
mon practice to write these equations as : 



k 2 F(k;0 



1 — — / rfk- 

np 2 J k 2 +M 2 (k;£) 

k 2 



1 



P 



2kp 



In 



k 



P 



k — p 



M(p;Q = a(g + 2) 

f(p;Q np 



kF^OMik^ 

k 2 +M 2 (k;0 



k + p 



where a — e 2 /An. Owing to the fact that in the Landau 
gauge, F(p; 0) = 1, it has long served as a favourite gauge 
where one only has to solve the following equation : 



M( P ;0) = — 
irp 



dk 



kM(k;0) 
k 2 +M 2 (k;0) 



In 



(2) 



The dynamically generated non perturbative solution of 
this equation is depicted as a solid line in Fig. QJ. It 
corresponds to the appearance of masses through self in- 
teractions. In this article, we shall apply LKFT on this 
dynamically generated solution to find the result in an 
arbitrary covariant gauge. As we are interested in the 
analytical insight alone for the time being, we try to rep- 
resent the solution of the FP in the Landau gauge with 
some simple analytical function. We note that for small 
momenta p, the solution for the mass function is roughly 
a constant. The large-p behaviour of the mass function 
can be studied by transforming Eq. (J2J into a differential 
equation with the appropriate boundary conditions [T^ . 
Note that in the integrand we can use the approximation 
ln\(k+p)/(k-p)\ = (2p/k)8(k -p) + (2k/p)0(p- k), 
which is valid for k » p and also for k « p. Sub- 
stituting this expression into Eq. J3J), we can re- write it 
in differential form as (we have rescaled the integral so 
that the corresponding solution of the differential equa- 
tion can be expressed in units of e 4 .) : 



d_ 

dp 



idM(p; 0) 



dp 



T M(p;0) = 



(3) 



with boundary conditions p 3 dM(p;0)/dp — 0\ p= a and 
M(p; 0) = Ojp^oo- The solution of this equation is 



M(p;0) 



1t z p 



C1J2 



2 



c 2 Y 2 



9 

7T^p 



,(4) 



where J(x) and Y(x) are the Bessel functions of the first 
and second kind, respectively. The second boundary con- 
dition imposes C2 = 0. As J(l/ y/p) — > 1/p when p — > 00, 
we conclude that A4(p; 0) — » 1/p 2 in this limit. The same 
conclusion is reached from the exact result, Fig. (JJJ, by 
noticing that in the numerical data, p 2 M(p; 0) is practi- 
cally constant near the upper end of the momentum re- 
gion under consideration, namely, 10 3 . Therefore, we ap- 
proximate the numerical solution by the following func- 
tion : 



M(p;0) 



M m 2 
p 2 + m 2 . 



(5) 



The numerical result suggests us to choose Mq = 0.10406 
and mo = 0.1 in order that the infrared and ultravio- 
let behaviour of the mass function are correctly repro- 
duced. With this approximation, we find that the re- 
sulting Fourier transforms involved in LKFT (as we shall 
see later) are still not trivial to perform. Therefore, we 
need a further simplification. Owing to the fact that the 
mass function has the following infrared and ultraviolet 
features : 



• region I : p << 

Mj(p;0) = M , 

• region II : p >> m 

. . , n Mn m 2 , 

Mi(p;0) = -V a 
p A 



(6) 



(7) 
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Mass Function in Landau Gauge in Rainbow Approximation 

Numerical Solution - 
Approximation 1 
Approximation 2 




FIG. 1: The mass function in the Landau gauge for the bare 
vertex. Approximations proposed in Eq. JSJ and Eq. JSJ are 
also shown. 



we can re-write the approximation to the mass function 
as follows : 



MO; 0) = M 



6(m - p) + ^§-9(p - mo) 



P 



(8) 



As shown in Fig. QJ, Eqs. 13181 provide a good approxi- 
mation to the mass function in the infrared and ultravio- 
let regions. We are now in a position to apply the LKFT 
to this approximate solution to find the corresponding 
result in an arbitrary covariant gauge. We carry out this 
exercise in the next section. 



III. LKFT AND THE DYNAMICALLY 
GENERATED MASS 

We start by putting forward the definitions and nota- 
tions we shall use along the way. We write out the FP 
in momentum and coordinate spaces, respectively, in its 
most general form as : 

S E (p;0 = Ato€) + t— r = y_ MM , (9) 

S E (x;0 = ^X(x;0 + Y(x;0 , (10) 

The subscript E stands for the Euclidean space. The 
above expressions are related through the following 
Fourier transformations 



Se(p;0 = d 3 xe ip - x S E (x;0 



Se(x;0 = J ^-e-^S E (p;0 



(11) 
(12) 



The LKFT relating the coordinate space FP in Landau 
gauge to one in an arbitrary covariant gauge reads 



Se{x;£) = S E (x;0)e~ 



(13) 



where a = a^/2. The way we proceed is as follows : We 
start with a dynamically generated solution of the SDE 
for the FP in the Landau gauge and Fourier transform 
it to coordinate space. We then apply the LKFT law. 
Fourier transform of this result back to the momentum 
space yields the FP in an arbitrary covariant gauge. 



A. Fermion Propagator in an Arbitrary Gauge in 
the Coordinate Space 

As we have seen in section II that in the Lan- 
dau gauge in the rainbow approximation, the func- 
tion F(p;0) = 1. Therefore, the FP can be written 
S E (p; 0) = l/iip 1 — M(p; 0)). Making use of this expres- 
sion in Eqs. (|1L)I12|> . the coordinate space FP in the Lan- 
dau gauge can be written as : 



X(x;0) 



Y(x;0) 



1 



2ir 2 x 2 



dp p 2 



1 

2^2 



p 2 +M 2 (p;0) 



dpp 2 M(p; 0) 
p 2 +M 2 (p;0) 



cospx 



sin px 



sin px 



px 



px 



(14) 
(15) 



For M(p; 0), we shall use the approximate solution given 
in Eq. JSJ. It naturally divides the integration region in 
Eqs. (|14I15|I into two parts. Therefore, 



X(x]0) = X l 
where 



X u and Y(x; 0) = Fj + Y u , (16) 



Xi 
X u 
Yi 
Yu 



1 



2ir 2 x 2 



dpp 2 
> 2 + M 2 



cospx 



smpx 
px 



1 




dp p e 




cos px - 


sin pa; 


2ir 2 x 2 




p 6 + Mlm\ 


px 




[•mo 


dp p 2 M 


sin pa; 








lo 


p 2 + M 2 


px 










dp p Mq TOq 


smpx 




"2^. 


I m 


p 6 + M§m 


4 



px 





(17) 
,(18) 
(19) 
(20) 



In order to be consistent with the approximations, 
Eqs. (I(il7l) made to arrive at Eq. JSJ), we neglect the corre- 
sponding terms in the denominators of the Eqs. H17I20JI . 
depending upon the region of integration over the mo- 
mentum space. Thus carrying out the radial integrations, 
we obtain 



Xi 

X u 
Yi 

Y n 



1 



2ir 2 M 2 x 5 



[3tooo;costoox + (? 



,2^,2 



3) sinmrja 



1 



47T 2 X 3 

mo 



7t + 2sinTOoa; — 2Si(rriox)] 
sin Too a; 



2ir 2 M x 2 
M m 
4tt 2 

Mo to 2 a; 
4tt 2 



cos moa; 



COS TOqX 



ITLqX 

sinTOoa; 



tuqx 



Si(mox) 



(21) 
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where Si(y) is the Sine Integral function, given by : 

P v sin t 

Si(y) = / V* ■ ( 22 ) 
Jo 1 



Now Eqs. (|1UI13I16I21|) constitute the FP in an arbitrary 
covariant gauge in the coordinate space. In the next sec- 
tion, we Fourier transform this result back to the mo- 
mentum space. 



B. Fermion Propagator in an Arbitrary Gauge in 
the Momentum Space 

In order to carry out the inverse Fourier transform, 
we use the auxiliary functions A(p; £) and £?(p;£) which 
are related to the more familiar functions M. (p; £) and 
F(p; £) in a simple fashion : 



M(p;0 = p 



2 MM) 



(23) 
(24) 



In terms of X(x;£) and F(x;£), the auxiliary functions 
can be written as follows : 



d 3 xe lpx Y(x;0 
4tt f°° 

— / dxx sin pxY(x; £) , 
V Jo 

d 3 xp- xe lp - x X(x;£) 



4n I dxx 2 
lo 



smpx 
px 



cospx 



(25) 



X(x;0- (26) 



We shall eventually be interested in the large and small 
limits of the momentum as compared to the dynamically 
generated mass scale to . However, as the coupling also 
has mass dimensions, special care may be needed for the 
evaluation of the the above integrals, especially for the 
S-term as we shall shortly see. We evaluate the A-teim 
and the i3-term separately below. 



1. The A-term 



Making use of Eqs. I|16i21jl . A(p;£) can be expanded 
term by term as : 



A( P ;0 

where 
Ai(p;0 = 



2m M m 

Ai{p; A 2 (p; £) 



ttMqP 
M m 2 
irp 



Tip 



A 3 (p;0 , 



(27) 



— e smpx 
o x 



cos rriQX 



sin niQX 
mox 



A 3 ( P ;0 



dxxe 



smpx 



cos toqX 



sin mox 



toqx 



dxx 2 e ax smpx 



- - Si{m x)\ . (28) 



The results of these integrations have been presented in 
the appendix. The addition and simplification of these 
expressions brings about some cancellations to yield the 
final simple result as : 



A(p;£) 



1 

7T 

a 
2?rp 



1 M m 2 (3a 2 -p 2 ) 
Mo {a 2 + p 2 ) 3 

1 M ml{a 2 -3P 2 ) 1 

w ~ 



\T+ + T_ 



(a 2 +p 2 ) 



2 v'i 



M m 
(a 2 + p 2 ) 2 



4a (3a 2 - p 2 ) 

h to — J - — g— 

7r a z +p^ 



(29) 



where 



T± — arctan 



to ±p 



L = In 



a 2 + (m +p) 2 



a 2 + (toq - p) 2 



(30) 



2. The B-term 



We carry out a similar procedure to evaluate B(p,£). 
However, the integrations involved are harder as we shall 
shortly see. Eqs. (|16I21(I permit us to write 



B(P,0 = ~^B 1 (p;0 + -\l 
+ B 3 (p;0--B 4 ( P ;0 , 



M 2 



(31) 



where the integrals involved are outlined below 



COS TOqX 



sm ttiqX 
nigx 



cospx - 



smpx 



px 



B 2 (p;0 
B 3 (p;£) 
B 4 ( P ;0 



dx 

x 
dx 

x 

dx 
— f 

x 



e ax sin toq x 



cospx 
cospx 



cospx - 

sinpx 

px 
sinpx 

px 



smpx 



px 



Si(mox) . (32) 



Except the last one, the other integrals are less compli- 
cated to evaluate. Making use of the same notations as 
before, these are listed below : 



t 



6tooP 
a 



24toqp 



[(to 3 +p 3 )T + -(to 3 ~p 3 )T_] 
[4to p + (a 2 + 3to 2 + 3p 3 )L] 



5 



mo 
2p 



[T- - T+] 



Ap 



-L 



ci p 
Bs(v;0 = — 1 H — arctan — 
p a 



(33) 



A word of caution is appropriate here. In order to verify 
our findings, we shall always want to make a connection 
of our results for the large and small values of momentum 
to the ones in the Landau gauge. It appears that B 3 may 
offer a tricky problem in this context because after taking 
the small p limit, we cannot retreat to the Landau gauge 
which corresponds to the a — > limit. The term B±(p\ £) 
comes to our rescue and it becomes evident if treated 
carefully. It is for this reason that we prefer to write it 
in the form : 



where 

In = 



JO 



dxx" 



n=0 



cospx 



sin px 
px 



(34) 



Si{niQx) . (35) 



This expression is perhaps the most convenient to deduce 
the FP for the asymptotic limits of momenta and then 
also to be able to make connection with the results of the 
Landau gauge. One can find out that when mo < p, 



In = 



m 
rtl+n 



r(l + n) p F q 
+ r(n) p F q 
and when mo > p 



1 n+l 



2' 2 
1 n+l n 



,1 + 77 



3 3 
2' 2 



2 ' 2 
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,(36) 



7T /7l7T\ r , „ , r(n) /717TN 

J ™ = ^ 008 (t) ^ n - ^ + r(n)1 + ^ 8,11 (t) 



C 1 
P r 9 



n + l n n 



1 

2' 



pFq 



n+lnnl f3 n I p' 
~ 2~' 2' 2 | ' 1 2 ,1 + ^ ' 5 ^ 



(37) 



Finally, Eqs. I|9I29I37() give the dynamically generated 
FP in an arbitrary covariant gauge. 



In order to make a quick comparison possible, we write 
below the expressions for A(p; 0) and B(p; 0) : 



1 



A ^ - "Mo 



0(m o - p) 



B(p; 0) 



P M 



p) 



9(p 



9(p - mo) 



mo) 



(38) 
(39) 



We now proceed to evaluate the asymptotic limits of the 
LKF transformed functions. 



A. Large-p Behaviour of the Fermion Propagator 

In the large-p limit, one can verify that 

Aam Q {ml + 3M§)~ 



-m Q M 



3M TT 



1 



B( P ;0 = -1 + 



(40) 
(41) 



The Landau gauge results, Eqs. (|38I39|) . are readily re- 
produced if we substitute a = 0. Correspondingly, the 
mass function and the wavefunction renormalization ac- 
quire the following form : 

4amo (ml + 3M 2 ) 



M(p;0 



mo Mo 



3M TT 



1 

p 2 



O 



1 



p- 



_ c 3 (0 



° ^p 3 



F(p;0 = 1 



■oil 



(42) 
(43) 



Again, obviously, if we set a = 0, we recover the Lan- 
dau gauge results. An important implication of these 
last equations is that in an arbitrary covariant gauge, 
the mass function will continue to fall as 1/p 2 for its 
asymptotically large values. Numerically, such an effect 
was seen although for a very small range of values of 
the gauge parameter (£ = — 1) in For a much 
broader range (£ = — 5), the same observation was re- 
ported in [lOj . Here, the analytical incorporation of the 
LKFT allows us to generalise these findings to an arbi- 
trarily large value of the gauge parameter. Moreover, 
the wavefunction renormalization continues to approach 
1 in the large p limit for arbitrarily large values of the 
gauge parameter. This fact is also in accordance with 
the numerical computations done for a limited range of 



IV. RELEVANT LIMITS 

We are now in a position to carry out the analytic 
analysis of the results obtained in the previous section. 
We analyse dynamically generated FP in the asymptotic 
limits of momenta, i.e., when p » mo and mo >> p. 



B. Low-p Behaviour of the Fermion Propagator 

An analogous analysis for the low-p regime yields 

a 4 V 7T 1 



6 



2arriQ 



•mg) 



7r(a 2 4 
_2 /J_ _ 

d(o + o(p 2 ), 



i 

Mo 



3to 2 M, 



arctan 
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m 



(44) 



3M 2 7T 



am (3a 2 + 5mg - 2M§) 



(a 2 + m„) 



2\2 



-3 arctan 



rn 



P 2 + G(p 3 



c 2 (Op 2 + o( p 3 ) , 



(45) 



where a neat cancellation has taken place between B3 
and the first sum in B4 (refer to the last equation of the 
Appendix). Taking the limit a — > 0, we again recuper- 
ate the corresponding results in Eqs. I|38I39|) . The mass 
function and the wavefunction renormalization can now 
be calculated to be : 



Expectedly, 



cm 



0(p 2 ), (46) 



-C 2 (0p 2 + O( P 4 ) . (47) 



• The mass function is flat for small values of p, and it 
falls off as 1/p 2 for its large values. This behaviour 
was observed for a very limited number of values of 
£ in 0, Hcj after a complicated numerical exercise. 
Here we prove it to be the case analytically in an 
arbitrary covariant gauge. 

• The wavefunction renormalization is also a con- 
stant for small values of p. For the large values, 
it approaches 1. This fact is not a feature of F(p 2 ) 
just in the neighbourhood of the Landau gauge. It 
is true independent of our selection of £. 

Thus the p-dependence of the dynamically generated 
FP for the small and large values of the momentum con- 
tinues to have the same qualitative features in an arbi- 
trary covariant gauge as the ones in the Landau gauge 
and in its neighbourhood, |jj . 



V. FULL VERTEX 

In the previous sections, we have relied on the assump- 
tion that the bare vertex is the correct approximation in 
the Landau gauge. Several works indicate that this is not 
the case. The latest in a series of proposals for the full 
vertex in QED3 is the one suggested in [||. However, its 
employment to solve the SDE for the FP is a formidable 
task even in the Landau gauge. In perturbation theory, 
as indicated in ^| (see Eq. (26) of that article), the 
transverse vertex in fact vanishes in the Landau gauge 




FIG. 2: Wavefunction Renormalization in the Landau gauge 
for the full vertex. Approximation proposed in Eq. 1481 is 
also shown. 



in the limit k 2 >> p 2 for QED in arbitrary dimensions. 
Based upon this, there are several suggestions of using a 
full vertex whose transverse part vanishes in the Landau 
gauge. In this article, we concentrate on all such ver- 
tices, e.g., |3 13 13 13- I n tms case > the numerical 
behaviour of the wavefunction renormalization modifies 
to the one shown in Fig. @. It is no longer 1 for all 
momentum regimes. Instead, it behaves like a constant 
(different from unity) for low momentum, and tends to 
one as p — > 00. Therefore, in addition to using the ap- 
proximation in Eq. © for the mass function (with Mq 
replaced by Mop and too by toof because the solutions 
in the Landau gauge for the bare vertex and the full 
vertex are not identical), we can use the following simple 
form for the wavefunction renormalization in the Landau 
gauge, Fig. (01 : 



F(p, 0) = F 0F 6{m F -p) + 0{p - m 0F ) . (48) 

Following the steps outlined earlier, it is easy to notice 
that only the terms X] and Yj get multiplied by the 
constant factor Fq, whereas, the rest remains unchanged. 
This fact makes the LKFT exercise straightforward. The 
slight modifications to the A(p;£ i ) and B(p;£) terms are 
the following : 



+ 
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2irp 



F, 
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M 0i ,m 2 F (3a 2 -p 2 )" 1 



M 0F 

Fqf 



(a 2 +p 2 ) 3 
M 0F ml F {a 2 



[T + 



3P 2 ) 



b(pA) 



M of 
M 0F m 0F 
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(49) 
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where B\, B 2 , -B3 and B4 have been already calculated. 
The large momentum limit for these expressions is now 



-o ^ 



A(p;0 = 



B(p;0 = -l + O - 
\P 



Consequently. 



AamoF(ml F F F + 3M 2 F ) 
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1 

¥ 

(51) 
(52) 



ml F M 0F + 



4am OF {ml F FoF + 3M 2 F ) 



SMqftt 



O 
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(53) 
(54) 



In a similar fashion, for the low-p behaviour, we now 
arrive at the following expressions for A(p; £) and B(p; £) 



A( P ;0 



m OF MoF 
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lam^F 
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( Fqf 
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7T V^/of a 4 
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arctan 



m 0F 



(55) 



B(p;0 = 



3M 2 f tt 



am 0F (3a 2 F F + 5m 2 F F 0F - 2M 2 F ) 



(a 2 + 771 2 ) 2 



-3-Fof arctan 



m 0F 



"OF 

■2 



P Z + 0{ P A ) 



c 2F (0p 2 + o( P 3 



(56) 



These functions easily get translated into the more famil- 
iar functions F(p;£) and Ai(p;(,) 



f(p;0 



c 2F (0 



(57) 

- C 2F {0P 2 + Oip*) . (58) 



We note that the incorporation of the full vertex does 
not change the qualitative features of the mass function 
and the wavefunction renormalization. The mass func- 
tion is flat in the infrared limit and falls off as 1/p 2 in 
the ultraviolet region. The wavefunction renormalization 
is also a constant for low values of p. It approaches 1 as 
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FIG. 3: Wavefunction renormalization for the bare vertex em- 
ploying LKFT. For a comparison, we also plot the results 
obtained by directly solving SDE. 



momentum p acquires large values. However, the con- 
stants multiplying these solutions get modified on em- 
ploying the full vertex. To get a quantitative estimate 
of this modification we carry out a numerical exercise in 
the next section. 



VI. NUMERICAL FINDINGS 

As stressed before, we cannot expect the numerical 
results presented in this section to be quantitatively ex- 
act and a true reflection of the choice of the full vertex 
owing to the facts that (i) the Landau gauge result has 
been approximated by a simple analytic function, and 
(ii) approximations I|tjl7f) have been employed in order 
to ensure that an analytical insight could be obtained. 
However, the approximate results indicate that we are 
making progress in the right direction. 

In order to view graphically our numerical results, we 
proceed as follows. From the corresponding expressions 
in the low and large momentum regimes for F and M, 
we perform the following parametrisation for the FP in 
arbitrary gauge : 



M{p;0 = M £{F ) 



6»(m c(F) -p) 



fle(p-m^ F} ) 



where 



F i{F )0(m i{F ) -p) + 0{p - m i{F) ) , 



(59) 
(60) 
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FIG. 4: Mass function for the bare vertex employing LKFT. 
For a comparison, we also plot the results obtained by directly 
solving SDE. 



and note that it obviously reduces to JSJ) and (|48l) on 

setting £ = 0. We separate the discussion on the bare 
and the full vertex in the following subsections. 



A. Bare Vertex 

In Fig. J3J) , we have plotted the wavefunction renor- 
malization in various gauges. Comparing them 
with the ones obtained by solving SDE, one sees 
that the difference is not enormous, reassuring the 
correctness of the method employed. 

In Fig. we have plotted the mass function in 
several gauges in order to compare the results of 
directly solving SDE against the ones obtained by 
employing the LKFT. With the increasing value 
of the gauge parameter, the direction of the shift 
in the mass function is diametrically opposed for 
both the cases. Which one should we trust more? 
LKFT arise as a necessary condition of gauge in- 
variance. Therefore, one would naturally regard 
the results more plausible after these transforma- 
tions have been taken into account. We shall see 
later that the studies carried out by including the 
full vertex justify this argument. 



B. Full Vertex 

SDE for the FP will yield an exact result if a full ver- 
tex is employed which complies with all the constraints 
imposed on it by the corresponding gauge field theory. 
Construction and usage of such a vertex is a prohibitively 
difficult task. Following the works of Ball and Chiu, we 
can divide the full vertex into two parts, the longitudinal, 
which ensures WGTI is satisfied and the transverse which 



FIG. 5: Wave function renormalization for the full vertex 
employing LKFT. 

Wavefunction Renormalization Comparison for the Full Vertex 



3. 



FIG. 6: A comparison of the results for the wavefunction 
renormalization (i) employing LKFT and (ii) results obtained 
by directly solving SDE for the full vertex. 



vanishes on contraction with the photon momentum vec- 
tor and remains undetermined by the WGTI. One of the 
most recent (and so far perhaps the most complete) at- 
tempts to propose such a vertex in quenched QED3 can 
be found in [6j. However, it is obviously a hard choice for 
its numerical implementation. For the purpose of being 
able to compare our results with earlier works we restrict 
ourselves to the family of transverse vertices which van- 
ish in the Landau gauge although an extension to a more 
general case is straightforward. We observe the following: 

• In Fig. (JSJ) , we have plotted the wavefunction renor- 
malization in several gauges and in Fig. ©, we 
present a comparison against the results of directly 
solving SDEs. The results are found to be in fairly 
good agreement. 

• In Fig. (JJJ, we have plotted the mass function in 
several gauges. In case of the full vertex, it is a bit 
hard to compare the results of directly solving SDE 



9 



LKFT-Mass Function in Various Gauges for the Full Vertex 



VII. CONCLUSIONS 



S o., - 



FIG. 7: Mass function for the full vertex employing LKFT. 
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The bulk of strong interaction phenomena require a 
non-perturbative approach. An example is the dynam- 
ical generation of fundamental fermion masses, which 
lies beyond the realm of perturbation theory. This phe- 
nomenon is linked with the non-perturbative behaviour 
of the FP and is governed by its SDE. As the SDE of 
the higher point Green functions are intricately related 
to it, a meaningful truncation scheme which maintains 
the key features of the theory, such as gauge invariance, 
continues to be a challenging problem. Along with other 
guiding principles, a correct inclusion of the LKFT of 
the FP and the fermon-boson interaction is crucial for 
arriving at reliable conclusions. We carry out a prelimi- 
nary study in this connection. Starting from a dynami- 
cally generated chirally asymmetric FP in one gauge, we 
perform an LKFT to find its form and properties in an 
arbitrary covariant gauge £. Encouragingly, we find ex- 
pected behaviour of the FP for the asymptotic regimes of 
momenta for all £, a fact supported by earlier direct nu- 
merical solutions of the SDE (where no reference is made 
to the LKFT) in a small region of £ close to the Landau 
gauge • The initial results for quenched QED3 pre- 

sented here provide the starting point for a more rigorous 
and exact numerical study. One can then go on to take 
up the more interesting case of QED4. The solution to 
this problem will in turn be the staring point for a study 
of QCD where the non-abelian nature of interactions, so 
essential for both confinement and asymptotic freedom, 
will further complicate the problem. All this is for the 
future. 



FIG. 8: A comparison of the results for the mass function (i) 
employing LKFT and (ii) results obtained by directly solving 
SDE for the full vertex. 
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against the ones obtained by employing the LKFT 
as the results in the former case are known only in a 
small region near the Landau gauge [3, Hfl uM ■ We 
offer a comparison in Fig. (8). An interesting thing 
we note is that in both the cases the trend for the 
function is to get shifted upwards with increasing 
values of the gauge parameter. This is the same 
behaviour as we had noted for the bare vertex by 
using the LKFT. Moreover, the difference between 
the two types of approaches seems to close down 
on each other when the full vertex is employed. It 
is understandable as the constraints of gauge in- 
variance are now being used in both the cases. A 
considerable advantage of using LKFT method is 
that the results for an arbitrary value of the covari- 
ant gauge parameter are readily available. 



Appendix 

We list below several of the integrals used, in this 
article for a quick reference. Fourier transforming the 
Landau gauge FP in the momentum space to the one in 
the coordinate space requires the following integrals : 



dpp 
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X 

2m 

x 2 



sin mx 



cos mx 



Si(mx) 



The Fourier trasnsform back to the momentum space af- 
ter having performed the LKFT makes use of the follow- 
ing results : 
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where 



In arrivng at the above relations, we have made use of 
the following relations 
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We have also made use of the following sum 
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